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ABSTRACT. This paper presents radar-altimeter scattering models for each of 
the diagcnctie zones of the Greenland ice sheet. A AFE radar-altimeter waveforms 
obtained during the 1991 and 1993 NASA multi-sensor airborne altimetry 
experiments over Greenland reveal that the Ku-band return pulse changes sig- 
nificantly with the different diagenctic zones. These changes are due to varying 
amounts of surface and volume scattering in the return waveform. 

Tn the ablation and soaked zones, where surface scattering dominates the A AFE 
return, geophysical parameters such as rms surface height and rms surface slope are 
obtained by fitting the waveforms to a surface-scattering model. Waveforms from the 
percolation zone show that sub-surface ice features have a much more significant effect 
on the return pulse than the surrounding snowpack. Model percolation waveforms, 
created using a combined surface- and volume-scattering model and an ice-feature 
distribution obtained during the 1993 field season, agree well with actual A AFE 
waveforms taken in the same time period. Using a combined surface- and volume- 
scattering model for the dry-snow-zone return waveforms, the rms surface height and 
slope and the attenuation coefficient of the snowpack are obtained. These scattering 
models not only allow geophysical parameters of the ice sheet to be measured but also 
help in the understanding of satellite radar-altimeter data. 


1 . INTRODUCTION 

It has been shown by several authors that the scattering 
and emission processes in the microwave frequency region 
are affected by sub-surface inhomogeneities in the 
Greenland ice sheet fZwally, 1977; Zvvally and Glocr- 
son, 1977; Comiso and others, 1982; Swift and others, 
1985). These sub-surfacc features have been recorded by 
glaciologists and are one of several diagenctic facies, or 
physical and chemical characteristics, used to divide ice 
sheets into zones. Tn 1962, Benson (1962) determined that 
the Greenland ice sheet has ablation, soaked, percolation 
and dry-snow zones. Changes in Ku-band satellite radar- 
altimeter waveforms over the different diagenctic regions 
of Greenland have been noted and analyzed (Ridley and 
Partington, 1988; Partington and others, 1989; Davis and 
Moore, 1993). These satellite waveforms show that 
returns from some areas of the ice sheet have both a 
surface- and volume-scattering component. The volume 
or sub-surfacc component is due to penetration of the Ku- 
band pulse into the ice sheet. Determining the amount of 
the radar return caused by surface reflection as opposed 
to sub-surface scattering is, however, very difficult. 
Another issue that needs more thorough analysis is the 
response of the radar to ice features in the percolation 
zone. A more in-depth understanding of the effects of the 


diagenctic zones on the radar-altimeter waveform will not 
only help determine if there arc biases in altitude 
measurements due to penetration but will also allow 
mapping of the diagenctic boundaries. 

Tn 1991, NASA conducted a multi-sensor airborne 
altimetry experiment over the Greenland ire sheet. This 
experiment was the first in a series of biennial experiments 
aimed at measuring small-scale changes in the topog- 
raphy and mass balance of the ice sheet and at 
understanding the effects of penetration of a Ku-band 
radar-altimeter pulse into the ice sheet. The NASA 
Advanced Application Flight Experiment (AAFE) Ku- 
band radar altimeter, refurbished by the Microwave 
Remote Sensing Laboratory at the University of 
Massachusetts in 1990, participated in this experiment 
along with the NASA Airborne Oceanographic Lidar 
, s (AOIa, the NASA Airborne Terrain Laser Altimeter 
System (ATLAS; and the NASA Ka-band Surface 
Contour Radar (SCR). The experiment consisted of 10 
flights over the Greenland ice sheet between 22 August 
and 20 September 1991 aboard the NASA P-3 aircraft. 
Positional information provided by three Global Position- 
ing System (GPS) receivers aboard the NASA P-3 was 
used to control the autopilot of the aircraft. This system 
allowed accurate positioning under the European Remote 
Sensing satellite (ERS-1) trajectory and enabled accurate 
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repeat missions on all flight lines. The missions were flown 
over the ice sheet from Sondrestrom Air Force Rase (now 
known as Kangerlussuaq) on the west coast of Greenland 
at an average altitude of 400 m above the surface, and the 
flight lines covered an area between 65° and 75° N. 

All of the remote-sensing radar and laser instruments 
provided information about the aircraft altitude above 
the ice sheet. The GPS receivers, on the other hand, 
calculated the aircraft altitude above the ellipsoid, which 
is a reference around the surface of the Earth. The 
difference between these two measurements is the actual 
ice-sheet altitude above the ellipsoid. 

In 1994, the AAFE radar altimeter, AOL laser 
altimeter and ATLAS laser altimeter participated in the 
Greenland experiment again, along with the Coherent 
Antarctic Radar Depth Sounder (CARDS) built by the 
University of Kansas. Between 23 June and 9 July 1993, 
the XASA P-3 aircraft flew 1 I missions that covered an 
area between 61° and 74 N on the ice sheet. 

The 1991 and 1993 experimental flights over Green- 
land passed over the ablation, soaked, percolation and 
dry-snow diagenetic zones of the ice sheet. Since the 
airborne altimeter is unable to distinguish between the 
reflective characteristics of the ablation and soaked zones, 
these regions are referred to as one in this paper. As the six 
representative waveforms in Figure 1 show, the AAFE 


radar-altimeter returns received from the various zones 
differed significantly in shape. Each of the waveforms in 
Figure 1 is normalized in power with a single-range bin 
equivalent to 2.77 ns or 41.67 cm, and the latitude and 
longitude where each waveform was received, along with 
the ice altitude above the ellipsoid at that point, are 
noted. These variations in the waveform shape are due to 
differences in the scattering properties of the various zones 
of the ice sheet. Figure 2 shows the boundaries of the 
diagenetic zones in Greenland as well as the location on 
the ice sheet where each waveform was obtained. 
Waveforms 1 and 2, which have a large slope or sharp 
leading edge and a sharp trailing edge, come from the 
ablation and soaked zone, respectively. The trailing edges 
on waveform 3, which is from the percolation zone, has a 
smaller slope than waveforms 1 and 2, while the leading 
edges remain sharp, and waveform 4, also from the 
percolation zone, has a notable second peak. Waveforms 
5 and G, which were obtained near the summit in the dry- 
snow zone, still have a sharp leading edge but their 
trailing edge is long and gently sloping. These airborne- 
altimeter waveforms do not look like typical satellite- 
altimeter waveforms due to differences in beamwidth and 
footprint size. These differences arc explained in greater 
detail in the discussion and concluding remarks section ol 
this paper. 
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Fig. 1. Six altimeter return ivaveforms obtained from the Greenland ice sheet during a September 1991 mission. 
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Fig. 2. Simplified distribution of the diagenelic zones in 
Greenland according to Benson, and the location of the six 
altimeter waveforms of Figure l . 


The interaction of the AAFE altimeter 13.9 GHz pulse 
with the different diagenelic zones can be determined by 
calculating the dielectric properties of the snow and the 
depth of penetration in each region. Since snow is a 
heterogeneous mixture of air, ice and liquid water, all of 
which have known dielectric constants, the dielectric 
constant of snow, f b5 can be approximated. Measurements 
in the microwave region show that, at a given frequency, 
e* depends mainly on the density of the snow, p s , and the 
percentage of liquid-water content, m v . As the percentage 
of liquid- water content increases, both the real, e' s , and 
imaginary, c /r s , parts of the dielectric constant increase. 
An increase in density also causes an increase in the 
dielectric constant. This increase, however, mainly affects 
the real part of c*. Dry snow, which has a liquid-water 
content of zero, has an ef that depends only on the snow 
density, and an c'f that depends on the snow density and 
on the dielectric constant of pure ice, fj, at the 
measurement frequency (Ulaby and others, 1986). 
Models such as the Debye-like model and the Polder 
Van Santeen model can be used to calculate the dielectric 
constant of snow. For snow that has some water content, 
the Debye-like model approximates the dielectric con- 
stant, c ws , as 


e' ws = 1 + 1.83 p„ + 0.02m'° 15 + 00 ' 3m v 3 ' 


i + (///., r 


^ ws — 


0-073(///oK 

1 + U/fof 


1.31 
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where the effective relaxation frequency of wet snow, is 
9.07GHz (Hallikainen and others, 1986). The Polder 
Van Santeen model approximates the dielectric constant 
of dry snow, e f is, as 


/a* = (1 + 051^ s ) 3 

f'\ is = 3.275 asc" 


( p, iis) 2 (2e , ,i s + 1) 


(e'i + 2f', ls )(e'j + 2(e',i,,r) 
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(Ulaby and others, 1986). Both models agree well with 


measured results (Stiles and Ulaby, 1980; Ulaby and 
Stiles, 1980; Hallikainen and others, 1986). Once the 
dielectric constant of snow at 13.9 GHz has been 
estimated, the attenuation coefficient, o, and the depth 
of penetration, fi p , at this frequency can be calculated 
using 

2tt 

a = —\IMs/7.\ (3) 

and 
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(Ulaby and others, 1982;. When the depth of penetration 
is known, the type of scattering in the ablation, soaked, 
percolation and dry-snow zones can be determined and 
used to model the altimeter return waveform. 


2. SCATTERING MODELS OF THE DIAGENETIC 
ZONES OF GREENLAND 


A. Ablation- and soaked-zones model 


Studies of the ablation and soaked zones in Greenland 
reveal that the upper layers of snow arc wet and dense 
due to summer melting. Measurements of snow density in 
the summer months range from 0.4 to 0.6 Mg m 3 
(Benson, 1962) and the percentage of liquid-water 
content may be as high as 3%. Estimating the dielectric 
constant of the snow using the Debye-like model and then 
calculating n and S p with Equations (3) and (4) gives an 
attenuation coefficient greater than 21Npm 1 or 
186 dBm 1 and a depth of penetration less than 2.0cm 
at 1 3.9 GHz. During the transitional season from the 
summer melting months of June and July to the winter 
season beginning around October, the ablation and 
soaked zones may not have a high water content but 
the snowpack will remain dense and still have a high 
attenuation coefficient. These approximate calculations of 
a and <*>,> show that there is little surface penetration at 
13.9 GHz, even in August and September, and therefore 
radar-altimeter returns in these regions are predomin- 
antly due to surface scattering. 

Surface scattering from a randomly distributed rough 
surface has been described by Brown (1977) as the 
convolution of the flat-surface impulse response with the 
distribution of height scalterers. The flat-surface impulse 
response, Pf s (f), is calculated by assuming a transmitted 
impulse and using the radar range equation 


pun = 


n A 2 r HI - f)GS0)aSf)) :l { 

(47t) J ilium arra /?■* 


(5) 


where the gain of ihe radar antenna is approximated as 

G{9) = G u c~ Al " 2(r ' /0! " (6) 

(Barrick, 1972). The backseat ter coefficient for the surface 
is 


o°{6) = 


HO*) - tan 2 0/s 2 

s 2 cos’ 6 


(7) 
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(Jackson and others, 1992), where T( 0°) is the power- 
reflection coefTieicnt at nadir and s is the rms surface 
slope, but, for small incidence angles, it can be 
approximated by 

«•(»)= (8) 

Convolving Pk(t) with a Gaussian height distribution of 
scattcrcrs with a standard deviation of cr^ leads to a 
rough-surface impulse response of 




Range Bins 


where 


Co 


P t A 2 G 0 2 r( 0°) 


32tt 2 
2 II 

T =t , 

C 

__ 2\/2ah 

2H ( 81n2 1 

f = ( — y + — 

C \ Oyf S 2 


-1 


( 10 ) 

( 11 ) 

(12) 

(13) 


and erfc is the complementary error function (Barrick, 
1972; Brown, 1977). The total rough-surface system 
response is then calculated by convolving the rough- 
surface impulse response with the system point-target 
response, which can be approximated as a Gaussian pulse 
with a 3dB width of r p and a standard deviation of 
rr p = 0.425r p (Brown, 1977). The result is the same as 
Equation (9), except t v is now defined as 

tp = v^^jJW- ( 14 ) 


The rms surface height, cr\ l7 and rms surface slope, .s, 
are the two surface parameters that determine the shape 
of the return waveform. A change in <7j, affects the slope of 
the leading edge of the return waveform, while a change 
in s affects the slope of the trailing edge. To determine 
values of rr h , ,s and H or range, AAFE return waveforms 
are fitted using the least-mean-squared error (LMSE) 
method (Carnahan and others, 19C9; Press and others, 
1988) to a non-linear five-parameter model, which is 
Equation (9) plus a parameter, a, representing the noise 
floor of the waveform, or 


= « + 


_C0_ C^,/f-> 2 c -(2ir/^> f f £p 

wv u 



Altimeter waveforms from the ablation and soaked 
zones have the same shape as the model in Equation (15). 
When the surface-scattering model is fitted to waveform 2 
of Figure 1, as shown in Figure 3a, the LMSE fit yields a 
a u of 0.1 2m and an s of 5.8\ In the percolation zone, 
however, return waveforms, such as waveforms 3 and 4 of 
Figure 1, do not look like typical surface-scattering 
returns. Waveform 4, shown in Figure 3b along with 
the model fit, has a trailing edge that differs considerably 
from a surface-scattering waveform, and therefore the 
values of s are not relevant. Although waveforms in the 
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Fig. 3. (a ) Surface-scattering model fit to waveform 2 of 
Figure I; (hj Surface-scattering model fit to waveform 4 
of Figure 1 ; (c) Surface-scattering model fit lo waveform 
6 of Figure /. 


dry-snow zone, such as waveform 6 in Figure 1, look like 
surface-scattering returns, filling them to the surface- 
scattering model, as shown in Figure 3c, results in 
abnormally high values of**. Il is unlikely that surface 
scattering is the only component contributing to the 
return waveform in the percolation and dry-snow zones. 

Even though the trailing edge of the waveform is 
sensitive to the different diagcnelic zones of the ice sheet, 
there is always a surface component in the return due to 
the dielectric interface between the air and the snowpack. 
This return from the surface creates the sharp leading 
edge shown in all of the waveforms in Figure 1. Since the 
rms surface height, <J\ l7 is determined from the leading 
edge of the return, it can be measured by fitting Equation 
(15) to AAFE waveforms. Figure 4a shows the ice 
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Fig. 4. Measurements of the rms surface height over a southwest to northeast flight line, (a) The actual altitudes measured 
and (b) the rms surface height determined by fitting 960 averaged waveforms to the surface-scattering model in Equation 
(15). 


altitudes from a southwest to northeast pass over the ice 
sheet and Figure 4b shows the corresponding rms surface- 
height values over the flight line, where each measure- 
ment represents a fit to 960 averaged waveforms. The 
average rms surface hei gin over this section of the ice 
sheet is 42 cm and the standard deviation from the mean ., 
is 33 cm. Figure 4b also shows the approximate location of 
the diagenelic transitions and the mean and standard 
deviation of the rms surface height for each individual 
diagenelic region of the ice sheet. The higher values for a\ i 
in the ablation and soaked zones correspond to the 
greater variability in ice altitude in Figure 4a, while the 
lower values in the dry-snow zone correspond to the 
smoother ice altitudes. 

Since the trailing edge of the altimeter waveform is 
affected by sub-surface and volume scattering, the rms~ 
surface slope, ,s, cannot be determined by fitting reluril 
waveforms from the percolation and dry-snow zones to a 
surface-scattering model. Returns from the ablation and 
soaked zones, however, are predominantly due to surface 
scattering and therefore can be used to measure s. Figure 
5a shows the ice altitudes over a section of the ablation 
and soaked zones and Figure 5b shows the corresponding 


rms surface-slope values resulting from fitting Equation 
(15) to AAFE return waveforms. Since a value of ,s 
greater than or equal to 15.6° represents a full antenna 
beam (see Discussion and Concluding Remarks for 
details), any .s greater than this is meaningless. Over 
tins flight section, the average and standard deviation of <s 
arc 9.0° and 3.6 C , respectively, but these slopes appear to 
be higher than expected. Since typical ocean rms slope 
values range between 5° and 8" (Jackson and others, 
1992), expected ice values for ,s would be even smaller. 
The trailing edge of the return waveform is also affected 
by aircraft mispointing angle and slope-induced error, 
which is the error induced on the return waveform due to 
a sloping surface within the radar-altimeter footprint 
(Brenner and others, 1983). Both mispointing angle and 
slope-induced error make the trailing edge of the 
waveform extend out farther and therefore create a 
larger value of s when the waveform is fitted to the 
surface-scattering model. Since the average mispointing 
angle over the flight section is less than 1.2 a and the more 
reasonable values of s correspond to the flat regions in 
Figure 5a while the high values correspond to the high 
slope areas, it is quite probable that the high values of .s 
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Fig. 5. Measurements of the rms surface slope over a flight 
section in the ablation and soaked zones, (a) The actual 
altitudes measured and (b) the rms surface slope 
determined by fitting 960 averaged waveforms to the 
surface-scattering model in Equation (15). 

arc the result of slope-induced error. For aircraft radar- 
altimeter measurements, these errors only affect the shape 
of the trailing edge and therefore, if they are not 
exceptionally large, they do not affect the range 
measurement, which is taken from the leading edge. 


B. Percolation-zone model 

In the percolation zone, where there is some summer 
melting, liquid water percolates down into the sub-surface 
snow and refreezes into ice layers, ice pipes and ice lenses. 
Stratigraphic studies reveal that these ice structures are 
located anywhere between 0.1 and 5 m below the icc- 
shect surface (Benson, 1962). Since this zone is a 
transition region between the soaked and dry-snow 
zones, the liquid- water content may vary from 3% near 
the border of the soaked zone to 0% near the dry-snow 
zone. Similarly, the density of the snow may range from 
0.5 to 0.3 Mg m \ Such a large change in the snow 
properties results in a vast range of attenuation co- 
efficients and depths of penetration in this region. 

The numerous odd-shaped sub-surface structures in 
the percolation zone, as well as the range of penetration 
depths, make it more difficult to model than the soaked 
and ablation zones. As the four altimeter return 
waveforms in Figure 6 show, the sub-surface ice features 
in the percolation zone affect the waveform shape. Since 
the ice lenses, pipes and layers are usually greater than 
10 cm in diameter and are large relative to the 2.16 cm 
altimeter wavelength, the resulting scatter from these ice 
structures is in the optical regime or the geometric optics 
limit (Ulaby and others, 1981). Therefore, the return 
from the ice features dominates any r scattering from the 
surrounding snow, and the percolation zone can be 
modeled as large volumc-scallcrers in a constant 
dielectric medium. As a result, the radar-altimeter return 
is calculated using the volume-scattering radar-range 
equation (Swift and others, 1985), which is 

„ P t A 2 r 2 r 6(t-2(f + ^))G 2 (6) 

— 3 / ft I 

(47T) J ilium vnl 11 

. c -2 „ A R-RJ vdV (16) 
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Fig. 6. Four sample altimeter return waveforms from the central part of the percolation zone of Greenland obtained during a 
September 1991 mission. 
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where the antenna gain G(0) is defined in Equation (6), T 
is the power-transmission coefficient, R 0 is the distance 
from the radar to the surface, R is a variable representing 
the range to the sub-surface scatterers and r s = c/fe^ is 
the velocity of propagation in the snow. Tins analysis is 
carried out in cylindrical coordinates and therefore dV = 
pdpd&dh where h is the depth below the surface. The 
extinction coefficient, k c , is defined as the sum of the 
power-absorption coefficient, and the power-scatter- 
ing coefficient, K s , or K e = k {1 + k 6 . Since the scattering 
losses in snow are negligible relative to the absorption 
losses at frequencies below 15 GHz (Stiles and Ulaby, 
1981; Hallikainen and others, 1986), k c for the percol- 
ation zone is 

Kc ~ = 2« (17) 

where a is the attenuation coefficient of the snow. Thus, 
the exponential component in Equation (16) represents the 
two-way power attenuation in the lossy dielectric media. 

The volume-backscatter coefficient, r] v> in Equation 
(16) is due to the ice lenses, pipes and layers and is defined 
as the product of an, the hackscatler cross-section per 
particle and n, the density per unit volume, or 

Vv = nan (18) 

(Swift and others, 1985). In the percolation zone, 
however, the ice features are not distributed evenly with 
depth and may be larger at certain depths. Therefore, 

n = n(p, (f>, h) = n(p, h) (19) 

ob = crii (p, 4>, h) = an(p, h) (20) 

assuming independence of (j). For a first-order approx- 
imation, the distribution is assumed to be independent of 
p resulting in n(h) and rrfj(ft). 

Substituting these results into Equation (16) and 
performing the integral gives 


where ,4 and B are the percentages of surface- and 
volume- ret urn power, respectively. 

To determine the volume-scattering component of the 
radar return, it is important to know the distribution of n 
and vs depth. Unfortunately, there is limited 
documentation on these parameters due to the difficulty 
in obtaining the measurements. Thus, in the 1993 
Greenland ground-truth experiment at Dye 2, a sample 
distribution of ice layers, lenses and pipes was obtained. A 
5 m by 3 m area of the ice sheet was probed at 10cm 
increments using a marked steel rod which was pressed 
into the snov pack until it hit ice. The raw data from this 
experiment are shown in Figure 7. To validate this 
probing experiment, a pit was dug alongside the 5 m by 
3 m area and each 10 cm vertical wall was documented 
for size and location of ice features. 

As Figure 7 shows, there was a continuous ice layer at 
a depth of 1 m and the thickness of this layer fluctuated 
between 1 and 3 cm. Oilier sample juts dug several 
hundred meters away from this pit also had an ice layer at 
a depth of 1 m. This layer may have been the surface from 
a previous summer and it shows that some ice layers in the 
percolation zone may extend over very large areas. The 
results from this probing experiment were used to create 
the sample distribution 'of ?j(/i) shown in Figure 8. It is 
clear from this sample that the ice features do not have a 
standard distribution. In this case, they tended to be 
located between 37 and 42 cm, between 46 and 49cm, 
between 59 and 61 cm, between 65 and 70 cm and 
between 78 and 81 cm. The ice layer at 1 m is not 
included in this distribution. 

Determining the backscatter coefficient a\\ is as 
difficult as calculating n(li). Figure 9 is a picture taken 
during the 1993 Greenland experiment of an ice pipe 
leading down to an ice lens. Since the AAFE altimeter is 
nadir-viewing, these ice features would be viewed by the 
instrument from their top surface. In the optical regime, a 
target with gradual curving surfaces and an arbitrary 
shape has a backscatter coefficient of 


-8 In 2 


P rv (r) =C,c n 


Tf T 


■f< 


r -8 In 2 c_ f 
T~ it 1 


" c - 2m >' 
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where C j and r are defined as 
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( 21 ) 

( 22 ) 

(23) 


In the percolation zone, the surface backscatter 
component of the altimeter return waveform is due. to 
the air snow 7 boundary, while the volume backscatter 
component is due to the optical volume scattering in the 
sub-surface snow. As a result, the total return power, 
P, (r), is the incoherent sum of these two components, 

Pr{r) = AP + PP rv (r) (24) 


an = 7T(1 1 a 2 1 r (0)| (25) 

where r(0) is the reflection coefficient at normal incidence 
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Fig. 7. A sample distribution of the ice features in a 5 m by 
3 m area of (he Greenland ice sheet at Dye 2. Each 
measurement represents the depth that the sleet rod hit ice. 
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Fig . 8 . A sample distribution of the location of ice features 
below the surface or n(h) in the Greenland snow pack at 
Dye 2 . 

and a i and arc the radii of curvature in the two planes 
of observation (Swift and others, 1985). Density measure- 
ments of the ice features taken during the 1993 
experiment were between 0.86 and 0.91 Mg m 3 . Since 
this is very close to the density of pure ice, which is 
0.9 16 Mg m 3 , the dielectric constant of pure ice, e*, may 



Fig. 9. Picture of an ice pipe and an ice lens taken during the 
1993 Greenland ground-truth experiment at Dye 2. 


be used to calculate r(0) or 


r(0) - 


y/fj- yfa 

yfei + 


( 26 ) 


Using ground-truth measurements of snow density and 
water content along with the Debye-like model in 
Equation (1) to approximate the dielectric constant of 
the surrounding snow results in e s values ranging from 
1.76 + j.001 to 1.83 4- j . 05 and an r(0) between 0.13 
j.006 and 0.17 j.008. 

Ice pipes in the percolation zone are typically 4 7 cm 
in diameter and their size is independent of depth in the 
first meter of the snowpack. In addition, since they are 
viewed by the radar altimeter from their top surface, they 
appear to be circular. Therefore, the average backscatter 
cross-section of an ice pipe in the percolation zone is 

ice pipe W — ice pipe ~ 28 dBsm. The ice lenses, on 

the other hand, are larger at greater depths. Lenses 
within 60 cm of the surface arc typically 10-20 cm in 
diameter, while those closer to 80 cm are often 50 60 cm 
in diameter. The resulting average backscatter cross- 
section for such ice lenses is (Tn] rv ] vnii {h = CO cm) « 
20dBsm and CTB\ vv \ niii (h = 80 cm) « — 8 dBsm. Finally, 
a continuous ice layer, such as the one throughout the 
V) e 2 area located at 1 m, is larger than the area of the 
AAFE footprint, and therefore its backscatter cross- 
section is limited by the footprint area to cr b layer 
(h — 1 m) « 20 dBsm. Since the ice feature with the 
largest contributing backscatter cross-section in the 
percolation zone is the continuous ice lay'er, it is the 
dominant seattercr, but the ice lenses and pipes are closer 
to the surface, and the incoming signal is not attenuated 
as much before it interacts with them. 

Although these ground-truth data give an idea of the 
size and distribution of ice features, they represent only 
one sample in the percolation zone, which covers a large 
area of the ice sheet. These results can be inserted into 
Equation (21), however, to get an idea of their effects on 
the return waveform. For example, a model optical 
volume-scattering return is created using Equation (21), 
the sample distribution of n{h) in Figure 8 along with the 
ice layer at 1 m and the ground-truth values of rr^. To get 
the total return power, P r , in Equation (24), this optical 
volume-scattering model is combined with a model 
surface return having typical values of 0.20 m and 4" for 
<j|, and .s, resulting in the final model percolation-zone 
waveform in Figure 10a. As expected, the ice layer has the 
largest effect on the waveform, creating the second peak 
at 1 m from the surface. Figure 10b is an actual AAFE 
return waveform obtained from the Dye 2 area during the 
1993 ground-truth experiment. This waveform exhibits a 
second peak similar to that in the model waveform in 
Figure 10a. The satellite altimeters do not typically 
observe these double-peak waveforms because the average 
depth of the ice layer varies over the much larger satellite 
footprint and therefore the ice-layer effect will not be as 
pronounced in the waveform. 

Both the model and actual waveforms agree with 
results from the ground-based Ku-band FMCW radar 
operated at Dye 2 during the 1993 ground-truth 
experiment. Nadir results for this radar show a sig- 
nificant return from the surface and from the ice layer at 
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Fig. 10. Altimeter model and actual returns from the 
percolation zone, (a) The combined surface and volume- 
scattering model and (b) actual AAFE return waveform 
from the percolation zone. 

I m and smaller returns from the ice lenses and pipes 
ijezek and others, 1994; Zabel and others, 1994; personal 
communication from I. Zabel, 1993). Therefore, when a 
continuous ice layer is present in the percolation zone, it is 
the dominant scatlerer and has the greatest effect on the 
A A F E re t u r n w a v e form . 

Although this model agrees with some of the altimeter 
returns from the percolation zone, it does not explain 
other more complex waveforms such as those in Figure 6. 
Since both n(h) and <7 b(/?) play an important role in the 
return waveform and both may vary significantly through 
the percolation zone, it is difFicult to invert the relation- 
ship in Equation (21) and determine n(h) and crn(h) from 
an actual waveform. Equation (21) docs, however, show 
the relationship between the return and n(h) and <T|;j (/?-), 
and allow sample distributions to be inserted and 
manipulated to illustrate the effect on the return 
waveform. 

C. Dry-snow-zone model 

In the dry-snow zone, where there is no summer melting, 
the density of the snow is low and the liquid-water 
content is 0%. Measured values of snow density, whichjie 
between 0.28 and 0.38 Mg m 3 (Benson, 1962) result in 
attenuation coefficients less than O.INpm 1 or 1 dBm’’ 1 
and depths of penetration greater than 5 m. Such a low 
attenuation coefficient suggests that a significant compon- 
ent of the return waveform is due to volume scattering 
from the sub-surface snow. 

The scattering due to the sub-surface volume, P rv (£), is 
again calculated using the volume-scattering radar-range 


equation (Swift and others, 1985): 

P,\ 2 T- f ^-2)(^ + ^G 2 W 

(47r) 3 /ilium vol 

. c -4o(R-R u )^ dF (27) 

The individual snow grains, which are less than 0.5 mm in 
this region, are small relative to the 2.16 cm wavelength of 
the altimeter. Therefore, the volume scattering is primarily 
due to Rayleigh scattering Ulaby and others, 1981). As a 
result, the volume-backscatter coefficient, ?/ v , is 

Vv = nan (28) 

where <Tb is the Rayleigh backseat ter cross-section per 
particle and n is the density per unit volume. In the dry- 
snow zone, the density of the snowpack and the size of the 
snow grains are fairly constant in the upper several meters 
(Ling, 1985), resulting in an a\\ and n that arc constant 
with position. The resulting return power due to the 
Rayleigh volume scattering is 

Prv (T) = c 2 - - 1 - - (c - 2 " r * T - e-*n r > 0 

pc — ZCVC-s 

p,. v (r) = 0 T < 0 (29) 

where 
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Equation (29) was obtained by performing the integral in 
Equation (27) in cylindrical coordinates. Davis and 
Moore (1993) obtained similar results by performing 
the analysis in spherical coordinates. Figure 1 1 shows how 
varying the attenuation coefficient, rv, of the sub-surface 
snow affects the shape of the volume-scattering return. As 
a approaches infinity, the depth of penetration goes to 
zero and, as expected, the volume-scattering return goes 
to zero. On the other hand, as r> goes to zero, the depth of 
penetration approaches infinity and P IV (r) becomes a 
function of the antenna beamwidth. 

As in the percolation zone, the total return power, 
P r (r), is the incoherent sum of the surface- and volume- 
backscatter components. The scries of waveforms in Figure 
12 illustrates the effect of varying the attenuation 
coefficient on the total return power. Tins is done by 
changing the snow density and liquid-water content using 
Equations (1) and (2). The first waveform shows that a 
liquid-water content of 3% results in a return that is 
predominantly due to surface-scattering. Sub-surface 
scattering effects do not appear until the liquid-water 
content approaches 0%. The following waveforms in 
Figure 12 illustrate the effect of keeping the liquid-water 
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Fig. 11 . Effects of varying attenuation , a , on the Rayleigh 
volume-scattering model of the dry-snow zone of Greenland. 


conlcnl al zero and reducing the snow density from 0.40 to 
0.27 Mgm 

Values of the attenuation coefficient in the dry-snow 
zone are determined using an LMSE fit to a non-linear 
seven-parameter model, which is the sum of Equation 
(24) and a noise-floor parameter or 

Pr( fit)(r) = a + APk{t) + BP n .{T). (34) 

This combined surface and Rayleigh volume-scattering 
model explains a return waveform such as waveform 6 in 
Figure 1, which has a sharp leading edge and a long, gently 
sloping trailing edge. As Figure 13 shows, the combined 
model fit to waveform 6 yields <7h = 0.42 m, s =2.6", o = 
0.1 Xpm 1 and <5 p =5.0m, which are reasonable surface 
and sub-surface parameters for the dry-snow zone. Using 




Fig . 12. Combined surface and Rayleigh volume-scattering 
content , m x . 
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Fig. 13. Combined surface and Rayleigh volume-scattering 
model fit to waveform 6 in Figure 1. 

the attenuation coefficient, a, and assuming the snow is dry 
and has m v =0%, results in an approximate snow density 
of Ps = 0.29 Mgm \ Attenuation coefficients from a section 
of the dry -snow zone obtained by fitting the combined 
model in Equation (34) to AAFE return waveforms have 
an average of 0 . 1 1 Xp m 1 or 1 dB m 1 . This is equivalent to 
an approximate snow density of = 0.28 Mgm Davis 
and Zw ally (1993) fitted a similar combined surface- and 
volume-scattering model to GEOSAT radar-altimeter 
waveforms and obtained values between 0.1 and 
0. 15 Xp m 1 in the dry-snow zone at 72 X, which was 
the upper extent of the SEASAT satellite coverage. 

Figure 1 4a and b show values of cr\i and s from the 

O 

same fit to Equation (34). The values oftrjj are similar to 
those shown in Figure 4b over this same flight section, 
which were obtained when the waveforms were fitted to 




lode! for varying values of snow density , p s and liquid-water 
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will create an AAFE waveform that is approximately 30 
40 range bins wide from the beginning of the leading edge 
to the end of the trailing edge. The SEASAT satellite 
radar altimeter had a similar pulse width of 3.2 ns, but its 
1.6° beam width and 800 km orbit (Davis, 1992) result in 
pulse-limited and beam-limited footprint radii of greater 
than 875 m and 1 1 km, respectively. As a result, SEASAT 
returns from the same rough surface will create wave- 
forms that are more than 100 range bins wide. The range- 
bin size of both systems is similar due to the equivalent 
pulse widths but this significant difference in the wave- 
form width or length of the trailing edge makes the two 
returns appear very different. 

Another cause for the differences between the airborne 
and satellite waveform shapes is the effects of rms surface 
slope, s, on the altimeter waveform. Figure 15a shows the 
effects of varying s on the AAFE return waveform for a 
constant altitude of 400 m. Once the rms surface slope 
approaches the same order of magnitude as the AAFE 
beamwidth, the waveform achieves its maximum width or 
the return is a full waveform. The occurrence of rms 
surface-slope values on the order of 15.6 n , however, is very 
unlikely and therefore surface scattering alone cannot 
produce a full waveform for the AAFE altimeter. As 
Figure 15b shows, the SEASAT returns approach a 


Fig. 14. Results of the combined surface and Rayleigh 
volume-scattering model fit to AAFE return waveforms 
from the dry-snow zone of the Greenland ice sheet . (a) 
Values of rms surface height , <r*, ; (b) Values of rms 
surface slope , .s. 

the surface-scattering model of Equation (15). This shows 
that, although Rayleigh volume scattering affects returns 
from the dry-snow zone, a surface-scattering model can he 
used for rms surface-height measurements, since the 
leading edge of the waveform is always due to a surface 
return. Values of .s, on the other hand, cannot be 
obtained using the surface-scattering model in the dry- 
snow zone. The combined model, however, allows part of 
the trailing edge to be attributed to s and part to a, 
resulting in reasonable values of rms surface slope, such as 
those shown in Figure 14b. This model, however, does not 
lake into account slope-induced error and mispointing 
angle of the aircraft, and therefore some of the data points 
may be biased slightly higher than expected. 

3. DISCUSSION AND CONCLUDING REMARKS 

As Figures 1 and 6 showed, the shape of the AAFE 
airborne-altimeter waveforms differs significantly from 
satellite returns from the Greenland ice sheet. One rca_son 
for this difference is the significantly smaller footprint of 
the airborne instrument. The AAFE altimeter has a 15.6° 
3dB beamwidth and a 2.77 ns compressed pulse width 
and, as a result, when the aircraft altitude is 400 m, its 
pulse-limited footprint is less than 20 m and its beam- 
limited footprint is less than 125 m. A reflective surface 
that is rough enough to reflect power back to the radar 
even at the farthest extent of the beam-limited footprint 




Fig. 15. Comparison of airborne- and satellite-altimeter 
waveforms. Effects of varying rms surface slope on (a) the 
airborne- A AF E-altimeter return waveform (assuming a 
constant altitude of 400 m and a beamwidth of 15.6°); 
(b) the satellite-altimeter waveform (assuming a constant 
altitude of 800km and a beamwidth of 1.6 °). 
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maximum width or a full waveform when the rms surface 
slope exceeds 1.6°, hut the AAFE returns will not be full 
until s exceeds 15.6°. Since the rms surface slope on the 
Greenland ice sheet is usually less than a few degrees, the 
effects of volume scattering on the trailing edge of the 
AAFE waveform arc easy to sec. The SEASAT satellite- 
altimeter returns, on the other hand, are typically full 
waveforms and the effects of volume scattering on the 
trailing edge are less evident. 

A third cause of the differences in waveform shape is 
that the area of the satellite-altimeter footprint is over 
2000 times greater than the AAFE airborne footprint. As 
a result, the deviation in surface altitude of the ice sheet 
across the satellite footprint can be much greater than 
across the airborne-altimeter footprint, causing a longer, 
more gently sloping trailing edge in the satellite-altimeter 
return. This also creates the double-ramp waveforms 
reported by Martin and others (1983), which were not 
recorded with the airborne AAFE altimeter. 

Despite the differences in the airborne and satellite 
waveforms, these models of the diagcnetic zones of 
Greenland provide a better understanding of how the 
radar-altimeter Ku-band pulse interacts with the ice 
sheet. In addition, these models allow important 
geophysical parameters to be measured, such as rms 
surface roughness and slope in the ablation and soaked 
zones, and attenuation and depth of penetration in the 
dry-snow zone. Although the percolation zone is complex 
and difficult to model, the ground-truth data taken at 
Dye 2 have helped in modeling this area. The next step in 
understanding scattering in this region is to use a 
multiple-scattering model to explain the complex radar- 
altimeter waveforms that are often obtained in the 
percolation zone. 

During the next few years, the NASA airborne group 
will repeat many of the original flight lines over Greenland. 
Comparisons between the 1991 and 1993 results and these 
future experiments will allow changes in the geophysical 
parameters of the ice sheet to be monitored. 
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